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Abstract 



^ ' An accurate study of the lowest lag and the low-lying excited lau, l'^u,g, ^^g,u electronic states 

of the exotic molecular ion -^3 in linear configuration parallel to a magnetic field is carried out. 
> , 
CO I The magnetic field ranges from 10^'^ G up to 4.414 x 10^^ G where non-relativistic considerations 

"^ ' are justified. The variational method is exploited and the same trial function is used for different 



magnetic fields. It is shown that the states of positive z-parity Icr^, Ivr^t, 1(5^ are bound states of 
the -^3"'' exotic ion for all magnetic fields studied. We also demonstrate that for magnetic fields 



Mh| B > 2.35 X 10^^ G the potential energy surface well corresponding to the lUg state contains at least 

o 

Vh ' one longitudinal vibrational state. It is also shown that the negative z-parity states lau, ^T^g, l^ui 

Ci ■ are purely repulsive in the whole range of magnetic fields studied, B = 10^'^ — 4.414 x 10^^ G. 
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I. INTRODUCTION 

Recently, it was announced that in a sufficiently strong magnetic field B > 10^^ G three 
protons situated on the magnetic line can be bound by one electron forming the exotic 
molecular ion H^^ in linear configuration 1]. A major feature of this configuration is that 
with a magnetic field growth the system becomes more and more bound (binding energy 
grows), and more and more compact (equilibrium distance decreases). Later it was shown 
that this exotic ion H^~^ in linear configuration becomes even the most stable one-electron 
system with smallest total energy for magnetic fields B > 10^^ G J2|. A recent study has 
demonstrated that the H^^ molecular ion can also exist in a certain spatial configuration - 
;he protons form an equilateral triangle while a magnetic field is directed perpendicular to it 
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The goal of this article is to make a detailed quantitative study of the ground state and 
the lowest excited states of the H^^-ion in a linear configuration for B = 10^*^— 4.414 x 10^^ G, 
which was announced in [Ij]. In particular, this is the ffist study of the lowest excited states 
of the H^'^ ion. 

Atomic units are used throughout {h='mf.=e=l) , although energies are expressed in Ry- 
dbergs (Ry). The magnetic field B is given in a.u. with Bq = 2.35 x 10^ G. 

II. GENERALITIES 

We study the system of three protons and one electron (pppe) placed in a uniform constant 
magnetic field. It is assumed that the protons are infinitely heavy (Born-Oppenheimer 
approximation of zero order) and that they are situated along the magnetic field direction 
forming a linear chain. The Hamiltonian which describes this system when the magnetic 
field is directed along the 2;— axis, B = (0, 0, B) is written as 

n=p' + ^ + ^+ I ------ + {pA + Ap)+A\ (1) 

i?+ K_ R+ + i?_ ri r2 r-^ 

(see Fig.^for the geometrical setting and notations), where p = —iW is the momentum, A is 
a vector potential which corresponds to the magnetic field B. It is chosen in the symmetric 
gauge, 

A = -^{-y, X, 0) . 

Hence, the total energy Et of the H^~^-ion is defined as the total electronic energy plus the 
Coulomb energy of proton repulsion. The binding energy is defined as the affinity of the 



system to form a bound state with respect to the system when the electron and the three 
protons are infinitely separated, E^ = B — Et- There are two dissociation processes: one 
of them has in the final state a hydrogen atom, Hf^ ^ H + 2p while the other one has a 
H2 molecular ion, H^^ -^ H2 + p- Therefore, the first dissociation energy is defined as the 
affinity of the system to form a bound state having two protons at infinity, E^^^^^ = E^ — Et , 
where Eh is the total energy of the hydrogen atom in a magnetic field B. While the second 
dissociation energy is defined as an affinity to form a bound state having one proton at 
infinity, E^^^^ = Ejj+ — Et , where E^+ is the total energy of the hydrogen molecular ion H^ 
in a magnetic field B. A contribution coming from the spin degrees of freedom is neglected. 
The problem we study is characterized by two integrals of motion: (i) the operator of 
the z-component of the angular momentum (projection of the angular momentum on the 
magnetic field direction) giving rise to the magnetic quantum number m and (ii) the spatial 
parity operator P{f -^ —r) which has eigenvalue p = ±1. Hence, any eigenstate has 
two explicit quantum numbers assigned: the magnetic quantum number m and the parity 
p. Therefore the space of eigenstates is split into subspaces (sectors) with each of them 
characterized by definite values of m and p. It is worth to note that the Hamiltonian (^J is 
also invariant with respect to z ^ —z. Therefore, in general one can classify the eigenstates 
using the quantum number cr = ±1 for positive/negative 2;-parity instead of p. However, 
this classification is related with the above described - there exists a relation between the 
quantum numbers corresponding to the z-parity and the spatial parity: 

p = a{-ir . 

To classify eigenstates we follow the convention widely accepted in quantum chemistry using 
the quantum numbers m and p. In particular, the notation for the states that we are going to 
use is similar to that introduced for H2-ion in parallel configuration |5] and is based on the 
following: the first number corresponds to the number of excitation - "principal quantum 
number", e.g. the number 1 is assigned to the ground state (lowest state), then Greek 
letters a, tt, 5 correspond to the states with m = 0, —1, —2, respectively, with subscript g/u 
(gerade/ungerade) describing positive/negative parity p. 

The excited states which we plan to study are the lowest states (of the type of the ground 
state) of the sectors with different magnetic quantum numbers m and p. It is quite obvious 
from a physical point of view that the ground state of a sector with m > always has larger 



total energy than those with m < 0. For this reason we restrict our consideration to the 
states with m = 0, —1, —2. 

Conservation of the z-component of the angular momentum implies that the electronic 
wave function (in cylindrical coordinates {p,ip,z)) can be written as 

^=e^-V™'^m, (2) 

where m is the magnetic quantum number. If we gauge rotate the Hamiltonian (^ with a 
factor from (J2I), it takes the form 

R+ R- R+ + R- ri r2 r^ 4 



where 



p^ = e-i^fp-M p ^im^p\m\ 



is the gauge rotated momentum (covariant momentum) and ipm are the eigenfunctions. The 
constant term mB describes the linear Zeeman effect splitting. In principle this term can be 
absorbed into the definition of the total energy. The representation is rather convenient 
since each Hamiltonian for fixed m describes the family of eigenstates with quantum number 
m and can be treated independently of the states with m' different from m. Now the 
Hamiltonian Q has only the invariance corresponding to the spatial parity conservation. 

As a method to explore the problem we use the variational procedure. The recipe of 
choice of trial functions is based on physical arguments and is described in full generality in 
p], where we address the reader. The ground state trial function for fixed m and p is chosen 
in the form 



with ipi . . .ipQ given by 



g-ai(ri+r2+r3)g-B/3ip2/4 if 0" = +1 



^i = { (5a) 

, if a = -1 

^2 = <^ ^ ^ (5b) 

0, ifcr = -l 

(^-a3{ri+r2) i g-QsCri+ra) i g-a3('-2+r3)'\ g-B/33P^/4 jf ^ _ i]^ 

V^3 = <( ' (5c) 

, if a = -1 

-a4r2-«5r3 , -«5r'2-a4r3'\ „-B/34P^/4 



_|_ ^g-"5»"2-a4r3\ 



e 



'g-a6{n+r2)-07'"3 _|_ Q-a(i(ri+r3)-a7r2 i g-a6{»'2+r3)-07ri'j g-B/^sP /4 j£ ^ _ i]^ 

i^5 = { ' ' ' (5e) 

) , if a = -1 

J. _ ('g-«8''l-"9''2-ai0''3 _|_ -^g-a9?'l-«8''2-Q:iO''3 _|_ g-ai0''l-"8''2-a9''3 _|_ ^C-p") 

-O8»'l-cnor2-a9»'3 _|_ „-cigr-i_-au)r2-as,r->, , ^g-«ior-i-«9r2-a8r3'\g-B/36P^/4 

where for the sake of convenience of a representation we use the quantum number a. The 
functions ipi . . .ipQ are 5*3 invariant with respect to the permutations of the identical protons 
and 0" = ±1 corresponds to the symmetric (antisymmetric) trial functions. It is clear that the 
functions ipi . . .i/jq are also eigenf unctions of the z-parity operator {pzipi = (Tipi, i = 1 . . . 6). 

In Eqs. (0), ^i,...,6 and ai^io,/?i,...,6, as well as the internuclear distances R± (see Fig. HJ 
are variational parameters [ij]. The total number of parameters is 23 for the symmetric 
case (o" = +1) and 10 for the antisymmetric case {a = —1) [1J|. The functions V'l.-.e carry a 
certain physical meaning. The function tpi {ip2) describes a coherent (incoherent) interaction 
of the electron with the protons. The functions ips, ip^ describe a coherent interaction of the 
electron with two protons (it may be called as H2 type interaction). In a certain sense it 
corresponds to the interaction oi H2 with a proton. The function ip^ represents an interaction 
of the electron with all three charged centers of the type {H2 — -ff-atom)-mixture. Finally 
the function ipQ is a nonlinear interpolation of the functions tpi . . .ipQ. 

Calculations were performed using the minimization package MINUIT from CERN-LIB. 
Numerical integrations were carried out with a relative accuracy of ~ 10^^^ by use of the 
adaptive DOIFCF routine from NAG-LIB. All calculations were performed on a dual PC with 



two Xeon processors of 2.8 GHz each. Every particular calculation of a given eigenstate at 
fixed magnetic field including minimization has taken several hours of CPU time. However, 
when the variational parameters are found it takes a few seconds of CPU time to calculate 
the variational energy. 

It is necessary to mention two technical difficulties we encountered. Calculation of two- 
dimensional integrals with high accuracy which appeared in the problem has required a 
development of a very sophisticated numerical technique. It was created a 'dynamical parti- 
tioning' of the domain of integration, which depends on the values of the variational param- 
eters (see e.g. j^). The domain partitioning was changed with a change of the parameters. 
Sometimes the number of sub-domains in particular integration was around 20. Another 
technical problem is related with a very complicated profile of the variational energy as 
function of the variational parameters which is characterized by many local minima, saddle 
points and valleys. Localization of the global minimum numerically of such a complicated 
function with high accuracy is a difficult technical problem which becomes even more diffi- 
cult in the case of twenty or more variational parameters. Examining the physical relevance 
of trial functions allows us to avoid spurious minima. The parameters obtained in (j3)) at 
every step of minimization were always examined from the physical point of view. Such 
considerations are always something of an art. 

III. RESULTS 

A. m = 

The m = subspace consists of two subspaces, p = 1 (even states) and p = —1 (odd 
states). 

1. lag state (p = 1) 

The state lag is the global ground state of the exotic H^'^-ion. Its eigenfunction has no 
nodes (Perron theorem). For this state our variational trial function (@)) with m = 0,p = 1 
depends effectively on twenty-two parameters (see footnote [13]). In comparison with the 
trial function used in the first studies jlll2] we have added two extra ansatz, ip4 and ip^. As it 
was mentioned above, the search for the global minimum numerically with high accuracy in 
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the case of so many variational parameters is a difficult technical task. We use a sophisticated 
strategy for localizing the minimum. As a ffist step we minimize ansatz by ansatz, then 
we take a superposition of two ansatz, then three ansatz etc. An essential element of the 
strategy is to impose a (natural) condition that the variational parameters change smoothly 
as a function of the magnetic field B. The above-mentioned strategy allowed us to improve 
our previous results for the ground state reported in [2| on total and binding energies, 
and also on the lowest vibrational energies. The qualitative results remain unchanged (see 
below) . 

The performed variational calculations indicate clearly to the existence of a minimum in 
the potential energy surface Et{R-\-,R-) of the (pppe) system for magnetic fields ranging 
B = 10^*^ — 4.414 X 10^'^ G. The minimum in the total energy always corresponds to the 
situation when R^ = i?_ = Req^ con&ming the qualitative symmetry arguments. Table 
m shows the results for the total Erp and binding energies Eh, as well as the internuclear 
equilibrium distance i?e/ for the ground state lag calculated with the trial function, 
and with the trial function used in our previous studies (Eq.(Q with A4 = A^ = 0) 
Thus, we predict that the exotic H^~^ ion can exist even in a larger domain of magnetic field 

riQ 

strengths as reported in [1|, |2|- We con&m a general statement that with an increase of the 
magnetic field strength the total energies increase, the system becomes more bound (binding 
energies increase) and more compact (the internuclear equilibrium distance decreases). 

In general, for all one-electron systems (in linear configuration) the binding energy in- 
creases asymptotically as oc log S for strong magnetic fields (for a discussion about the case 
of the hydrogen atom see e.g. p])- In the domain we study (B = 10^*^ — 4.414 x 10^^ G) the 
rate of increase of binding energy for H^'^ is slightly larger than the corresponding rates for 
the H atom and for the H2 ion (see j^|5|). As a result, for magnetic fields B > Ax 10^" G the 
H^~^ ion becomes more bound than the H atom, and for magnetic fields B >3x 10^^ G, H^'^ 
becomes the most bound (having the lowest total energy and the largest binding energy) 
among the one-electron systems H, H2, H^^ ■ 

An straightforward analysis of the internuclear equilibrium distances shows that in the 
domain B = 10^'^ — 4.414 x 10^^ G the rate of decrease of the internuclear equilibrium 
distance for the H^^ ion is also larger than the corresponding rate for the H2 ion (i?e/ 



decreases in ~ 20 times, while Req^ decreases in ~ 10 times when we go from B = 10 G to 
4.414 X 10^^ G). Making a comparison between i^l"*" and H2 one can see that the internuclear 



equilibrium distance for iJ|^ is always larger than that for H2- The internuclear distances 
converge to each other as the magnetic field increases: Req = 0.110 a.u., while Req = 0.102 
a.u. at 5 = 4.414 x lO^^ q (^ggg ^^^)_ 

Another important feature of the system is the behavior of its dissociation as a func- 
tion of the magnetic field. There are two dissociation processes: H^'^ — > H2 + p (i) and 
Hf^ ^ H + 2p (ii). Table El shows the dissociation energies corresponding to these processes 
Ed^^^ = E^+ — Et and -Edatom = ^h — Et , respectively, for different magnetic fields. Let us 
consider the first process (i). A negative dissociation energy Ed^^^ indicates that the ion H^'^ 
is unstable towards the decay H^^ -^ H2 + p- In particular, for B > 10^° G Ed^„^ is negative 
and decreases with a magnetic field growth reaching the minimum at i? ~ 10000 a.u. Then, 
for further magnetic field increase the dissociation energy (which is still negative) starts to 
increase monotonously. Eventually, for a magnetic field 5 ~ 3 x 10^^ G the dissociation 
energy becomes zero which means Et = E^+ and then starts to be positive. It implies that 
the ion H'^'^ becomes more bound than H2 and, in fact, this ion becomes the most stable 
one electron system. The above described behavior for the dissociation H^^ -^ H2 + p also 
indicates that, for a broad domain of magnetic fields, there exist two different values of the 
magnetic field for which we have the same dissociation energy Ed^^^- 

The situation is different for the second dissociation process (ii), Hl'^ ^ H + 2p. The 
dissociation energy Ed^^^^ increases monotonously in all the range of studied magnetic fields. 
In the domain 10^*^ — 4 x 10^° G the dissociation energy Ed^^^^ is negative and it indicates 
that the ion H^^ is unstable towards the decay H^^ -^ H + 2p. For magnetic fields 
i? > 4 X 10^° G the dissociation energy is positive and the ion H^'^ is more bound than the 
H atom. For comparison the dissociation energies -Edatom-ion fo^ the process H2 ^ H + p are 
also shown in Table IJ [l^. For magnetic fields B > 10^^ G, the dissociation energy -Edatom 
for the dissociation H^'^ — > if + 2p is smaller than the corresponding dissociation energy for 
H2 ^ H + p except for the domain B > 10000 a.u. 

The improvement in the total energy obtained using the trial function ffl in comparison 
to the results based on the reduced trial function (JD) {A4 = A^ = 0) jj, |5[ is the order of 
~ (1 — 5) X 10^'^ Ry for the whole range of magnetic fields studied. It represents a relative 
improvement of ~ 0.01% — 0.05% in the binding energies (see Table HJ. 

We should notice that in our previous studies |l|, I4I a different definition of the magnetic 
field atomic unit was used for unit conversion [Bq = 2.3505 x 10^ G). It leads to a relative 
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difference in the total energies of order of 10 ^ . It should be taken into account when a 
c„,npa.so. of the .esult. in na and the p..ese„t ...* . .ade. 

We show in Fig.|21the behavior of the variational parameters of the trial function (@)) as a 
function of the magnetic field strength. In general, the behavior of the parameters is rather 
smooth and very slowly-changing, even though the magnetic field changes by several orders 
of magnitude. In our opinion such a behavior of the parameters of our trial function Q 
reflects the level of adequacy (or, in other words, indicates the quality) of the trial function. 
In practice, the parameters can be approximated by the splines method and then can be 
used to study magnetic field strengths other than those presented here. 

a. Potential Energy Surfaces, Vibrational Energies. We carried out a detailed accurate 
study of the electronic potential surface Et{R+, R^) around the minimum and also along 
the valleys and around the barriers which ensure the existence of bound states. It allowed 
us to estimate accurately the (lowest) longitudinal vibrational energies. 

Let us first proceed to a description of the valleys of potential energy surfaces. Every 
potential surface is characterized by two valleys originated from the minimum. Those valleys 
are symmetric with respect to the bisectriz i?+ = i?_. Therefore it is sufficient to study a 
single valley and in further considerations we will be focused on one of the valleys which is 
almost horizontal. In Fig. |3] the valleys in the (i?+,i?_) plane and the profile along a valley 
are shown for different magnetic fields. Every profile is characterized by a minimum which 
corresponds to the equilibrium {R^ = i?_ = R^q) and a potential barrier. The height of 
the potential barrier is defined with respect to the energy corresponding to the equilibrium 
position (minimum) AEt = Eip""^ — ii^™". The asymptotics of the profile corresponds to 
decay of the system to H2 + p, for example, when i?+ -^ 00 (and -R_ — > Req^ ) or similarly 



R- —>■ 00 (and i?+ —> Re^ 



?^2 
•'eg 



The pattern of the valley in the (-R+, R-) plane exhibits a rather complicated behavior (see 
Fig. ^ where the behavior near minimum is displayed in amplified form). For example, for 
the horizontal valley the R- firstly grows reaching some maximum, then decreases reaching 
a minimum and after that increases approaching to asymptotics in i?+ from below. A similar 
picture holds for all magnetic fields studied. We do not think that this type of behavior is 
an artifact of insufficient accuracy of our calculations. In principle, the energy profile curve 
allows us to estimate the lifetime of the system with respect to the decay H^^ -^ H2 + p 
for magnetic fields 10^° G < i? < 4 x 10^^ G. However, such study is beyond the scope of the 



present analysis. 

Table HT] shows the values of the height of the potential barrier {AEt = E^"-^ — EJp^"-) 
along the valley of total energy and the position of the top of the energy barrier. From this 
table one can see that for the domain of magnetic fields B = 10^^ —4.414 x 10^^ G the height 
of the potential barrier increases dramatically from ~ 0.01 Ry up to ~ 5 Ry, respectivelv|l6l|. 

Now let us proceed to a calculation of the lowest longitudinal vibrational state. The 
method we use is based on the harmonic approximation of the potential around the min- 
imum. All necessary definitions to perform the analysis are similar to those which are 
used to study the vibrational states of linear triatomic molecules and can be found in 
standard textbooks (see, for example, l9|). Following the settings of FigHJ we find eas- 
ily the normal coordinates: Rg = 4=(i?+ -|- i?_) (it corresponds to the "symmetric" nor- 
mal mode), and Ra = -j^{R-^^ — R-) (it corresponds to the "anti-symmetric" normal mode). 
The lowest vibrational energy is then calculated as Eq'^'^ = Eq + Eq, where Eq = -i /^ and 
Eq = ^/f^ are the lowest "symmetric" and "anti-symmetric" vibrational energies respec- 
tively, rUp = 1836.15 is the proton mass measured in units of the electron mass, and kg, 
ka are the curvatures taken in a.u. near equilibrium for the symmetric and antisymmetric 
modes, respectively. 

The results for the lowest vibrational energies Eq^'' for different magnetic fields are shown 
in Table im These results indicate that with a magnetic field growth the form of the potential 
well around the minimum becomes sharper and the lowest vibrational energies increase 
drastically growing from ~ 0.1 Ry at 5 = 10^^ G up to ~ 2 Ry at S = 4.414 x 10^^ G PJ. 
Then, we can conclude that the energy well contains at least one vibrational state for 
magnetic fields B > 1000 a.u. 

b. Electronic distributions. The equilibrium configuration is characterized by a one- 
peak electronic distribution centered around the middle proton, being drastically shrinked 
in the direction transverse to the magnetic field as compared with the longitudinal direction 
(see Fig. it is worth to note that the different scales are used for the x and z axes). 
This one-peak form of the electronic distribution for equilibrium position is found for all 
magnetic fields studied. Also the electronic cloud has always a needle-like form. When 
the system is out of equilibrium, the electronic distributions show the appearance of a 
pronounced shoulder which follows the position of the more distant proton, reducing its 
height as the system approaches to the decaying configuration, where the electronic density 
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TABLE I: Total E^ and binding Ef, energies, internuclear equilibrium distances i?e„ and dissocia- 



tion energies E^,^^ = E^+ - Et (^3+ ^ H+ +p) and Ea^^^^ = Eh - Et {H^^ ^H + 2p) for the 



state 1(T„. Dissociation energies Ed^. 



Eh — E^+ {H2 -^ H + p) are shown for comparison. 



Results marked by asterisks (*) correspond to a trial function Q with ^4^5 
function, Turbiner and Lopez- Vieyra, unpublished 2003). 

B Et (Ry) E^ (Ry) 



(four Ansatz 



iJeg (a.u.) Ea,^^ (Ry) 



^d«o„ (Ry) 



Ed , ■ (Ry) 

"atom-ion V J / 



IQiOG 


1.8424 


2.4129 


2.072 




1.8438 


2.4116 


2.061 


10 a.u. 


6.6084 


3.3916 


1.431 




6.6094 


3.3906 


1.429 


10" G 


36.4297 


6.1234 


0.801 




36.4327 


6.1205 


0.799 


100 a.u. 


91.3611 


8.6389 


0.579 




91.3655 


8.6345 


0.578 


1012G 


410.3739 


15.1580 


0.345 




410.3762 


15.1558 


0.345 


1000 a.u. 


979.2171 


20.7829 


0.259 




979.2206 


20.7794 


0.259 


10" G 


4220.9286 


34.3905 


0.166 




4220.9320 


34.3872 


0.166 


10000 a.u. 


9954.5918 


45.4082 


0.130 




9954.5972 


45.4028 


0.130 


4.414 X IO13 G 


18727.7475 


55.2312 


0.110 




18727.7521 


55.2267 


0.110 



-0.7519 



-0.9581 



-1.3865 



-1.6521 



-1.9845 



-1.9956 



-1.3634 



-0.3890 



0.7294 



-0.2020 



-0.1039 



0.4105 



1.0596 



3.2427 



5.4615 



11.8395 



17.4532 



22.7305 



0.5500 


present 
(*) 


0.8542 


present 




(*) 


1.7966 


present 




n 


2.7117 


present 




n 


5.2272 


present 




n 


7.4571 


present 




n 


13.2029 


present 




n 


17.8422 


present 




n 


22.0010 


present 




(*) 



rearranges to mimic the electronic distribution of the H2 ion. Nothing especial seems to 
occur for the configuration near the top of the potential barrier. To illustrate the above- 
mentioned features we present on Fig. El the normalized electronic density distributions 
\E'^(x, y = 0, z)/ J \E'^(x, y, z)df and the corresponding contour distribution for the magnetic 
field B = 10000 a.u. for different positions of the system H^^ along the valley in potential 
surface in R^,R_. In this presentation we always keep the central proton at 2; = 0. The 
following points along the valley are chosen: (a) the equilibrium position, R^ = i?_ = 
0.13 a.u., (b) a point between positions of minimum and maximum, i?+ = 0.23 a.u., (c) a 
point close to the maximum, i?+ = 0.35 a.u. (maximum is around R+ ^ 0.36 a.u. - see 
Table |Tl)) and (d) a point on the asymptotic part of the valley where one of the protons is 
practically detached from the system. 
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TABLE II: Energies of the lowest longitudinal vibrational states (-E"''), height of the potential 
barrier (AEt = EJp""^ — E"™") and position of the maximum {R^(Emax) , R~ {Emax)) in the valley 
corresponding to the ground state lag of the ^3^ molecular ion in a magnetic field. Results marked 
by asterisks * are interpolated values of the height of the barrier using a cubic polynomial fit in the 
variable log(l + -B^) and the position of the maximum using a fit of the form R±{Emax,B) = I/P2 , 
where Pg is a quadratic polynomial in the variable log(l + B'^). 



B 






E^o'" (Ry) 


AEt (Ry) 


R^ 


-{Emax) (a.u.) 


R- 


-{Emax) (a.u.) 


lO^i G 






0.0951 


0.0101 




0.94 




0.758 


100 a.u. 






0.154 


0.1007* 




0.77* 




0.517* 


10^2 G 






0.338 


0.267 




0.58 




0.294 


1000 a.u. 






0.530 


0.581 * 




0.50* 




0.221 * 


10^3 G 






1.088 


1.895 




0.40 




0.147 


10000 a.u. 






1.594 


3.339* 




0.36* 




0.118* 


4.414 X 10^ 


^3G 




2.078 


4.815 




0.33 




0.101 


2. l(Ju 


state 


(P-- 


= -1; 













We have performed a detailed study for the state IcTu of the H^^ ion in symmetric 
configuration i?+ = i?„ = R in the domain of magnetic fields B = 10^'' — 4.414 x 10^^ G. 
The calculations were carried out using the trial function (@]) with m = and p = —1 
corresponding to a state of negative z-parity {a = —1). For a broad domain of internuclear 
distances R the electronic potential curve Ex = Et{R) shows no indication of the existence 
of a minimum. Hence, it is natural to conclude that this state is a purely repulsive state for 
B = 10^° — 4.414 X 10^^ G. It means that the molecular ion H^^ does not exist in the lau 
state as a bound state. 



B. m 



The subspace consists of two subspaces, p = 1 (even states) and p = —1 (odd states). 
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1. ItTu state (p = —1) 

We study the l7r„ state of (pppe) system in symmetric configuration Rj^ = i?_. = R. Tlie 
trial function (j3]) is taken with m = —1 and p = —1 (which imphes a = 1). The variational 
calculations indicate that the potential energy curve Et = Et{R) has clear minimum for 
magnetic fields B = 2.35 x 10^° G - 4.414 x 10^^ G. It manifests that if|+ ion exists in Ivr^ 
state. The results are presented in Table ITTTl Similar to the state lag the binding energy of 
the liTu state grows steadily with a magnetic field increase while the equilibrium distance 
shrinks in a quite drastic manner. For small magnetic fields the state Ivr^ is more extended 
than lag. The equilibrium distance Reg for ItTu is much larger than this distance for the lag 
state, as for large magnetic fields these equilibrium distances become comparable. Among 
the m = — 1 states the state iTr^ has the smallest total energy. A comparison with the 
corresponding energy of the Ivr^ state of H2 (cf. Table VII in [5]) shows that the energy of 
the l7r„ state of H^^^ is always larger for B = 2.35 x 10^° G - 4.414 x 10^^ G. 

Transition energy from the lag state to the Itt^ state, Eit,^ — Ei^g is easily calculated 
by taking data from Tables U and IIIIl In the whole range of magnetic fields studied the 
transition energy increases monotonically with the magnetic field growth as expected (see 
Fig. ED. 

In Fig. Efa) a plot of the normalized electronic density distribution 
\E'^(x, y = 0, z)l j \l'^(x, y, z)dr for the ItTu excited state in a magnetic field oi B = 10000 a.u. 
is shown. By looking at the corresponding contour distribution it is evident that the elec- 
tronic cloud has a torus-like axially-symmetric form with respect to the 2;-axis. The radial 
size of the torusSlS] decreases with a magnetic field increase. A similar qualitative behavior 
of the electronic density distribution is observed for different magnetic fields . 

2. 1-Kg state (p = 1) 

We have performed a detailed study for the state Ivr^ of the H^'^ ion in symmetric 
configuration i?+ = i?_ = R in the domain of magnetic fields B = lO^*' — 4.414 x 10^^ G. 
The calculations were carried out using the trial function (0)) with m = —1 and p = 1 
corresponding to a state of negative z-parity {a = —1). For a broad domain of internuclear 
distances R the electronic potential curve Et = Et{R) shows no indication of the existence 
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TABLE III: Total Ex, binding Ei, energies and equilibrium distance Reg for the excited state l7r„. 

B^ Et (Ry) Eb (Ry) Reg (a.u.) 

lOa.u. 7.9289 2.0711 2.413 

10^1 G 38.6589 3.8943 1.238 

100 a.u. 94.3927 5.6073 0.869 

10^2 G 415.3661 10.1658 0.497 

1000 a.u. 985.7952 14.2048 0.366 

10^3 G 4231.0542 24.2649 0.226 

10000 a.u. 9967.3669 32.6331 0.174 

4.414 X 10^3 G 18742.7564 40.2223 0.145 

of a minimum. Hence, it is natural to conclude that this state is a purely repulsive state for 
B = 10^'' — 4.414 X 10^^ G. It means that the molecular ion H^^ does not exist in the Ivr^ 
state as a bound state. 
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C. m = -2 

The subspace consists of two subspaces, p = 1 (even states) and p = —1 (odd states). 

1. 15g state (p = 1) 

We study the 16g state of (pppe) system in symmetric configuration i?+ = _R_ = R. The 
trial function Q is taken with m = —2 and p = 1 (which imphes a = 1). The variational 
calculations indicate that the potential energy curve Et = Et{R) has clear minimum for 
magnetic fields B = 2.35 x 10^° G - 4.414 x 10^^ G. It manifests that if|+ ion exists in l^^ 
state. The results are presented in Table HVl Similar to the states lag and Ivr^ the binding 
energy of the 16g state grows steadily with a magnetic field increase while the equilibrium 
distance shrinks in a quite drastic manner. For small magnetic fields the equilibrium distance 
R = Req for 15g is much larger than this distance for the lag and Itt^ state, as for large 
magnetic fields all three equilibrium distances become comparable. Among the m = —2 
states the state 16g has the smallest total energy. A comparison with the corresponding 
energy of the ISg state of H2 (cf. Table IX in [5|) shows that the energy of the 15g state of 
H^+ is always larger for B = 2.35 x 10^° G - 4.414 x lO^^ G. 

Transition energies from the ItCu state to the 16g state Eis —Ei^^^ can be easily calculated 
from Tables ITTTl and HVl showing a monotonically increasing behavior in all range of magnetic 
fields studied (see Fig. E)). It is worth to note that the transition energy from Itt^ state to 
the ISg state is always smaller than the transition energy from lag state to the Ivr^ state. 

In Fig. Efb) a plot of the normalized electronic density distribution 
\l/^(x, y = 0, z)/ f \l/^(x, y, z)dr for the 15g excited state in a magnetic field oi B = 10000 a.u. 
is shown. Similar to the case of the Ivr^ state, the electronic cloud for the 15g state has 
a torus-like axially symmetric form with respect to the magnetic field line (z-axis). It is 
worth to note that the radial size of this torus-like form of the electronic distribution for 
the 15g state is larger than for the l7r„ state. 

2. 15u state (p= -1) 

We have performed a detailed study for the state iSu of the H^'^ ion in symmetric con- 
figuration i?+ = R_ = R in the domain of magnetic fields B = 10^° — 4.414 x lO^'^ G. 
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TABLE IV: Total Et, binding Ef, energies and equilibrium distance R^q for the state 15g. 

_B Et (Ry) Eb (Ry) Reg (a.u.) 

lOa.u. 8.3644 1.6356 3.214 

10^^ G 39.4450 3.1082 1.548 

100 a.u. 95.4919 4.5081 1.070 

10^2 Q 417.2496 8.2823 0.601 

1000 a.u. 988.3293 11.6707 0.437 

10^3 G 4235.0826 20.2366 0.266 

10000 a.u. 9972.5387 27.4613 0.202 

4.414 X 10^3 G 18748.9067 34.0720 0.167 

The calculations were carried out using the trial function Q with m = —2 and p = —1 
corresponding to a state of negative z-parity {a = —1). For a broad domain of internuclear 
distances R the electronic potential curve Et = Et{R) shows no indication to the exis- 
tence of a minimum. Hence, one can conclude that this state is a purely repulsive state for 
B = 10^'^ — 4.414 X 10^^ G. It means that the molecular ion H^^ does not exist in the ISu 
state as a bound state. 



IV. CONCLUSION 

We have carried out an accurate, non-relativistic calculation in the Born-Oppenheimer 
approximation for the low-lying states of the H^'^ molecular ion in the linear parallel con- 
figuration placed in a constant uniform magnetic fields ranging from B = 10^*^ G to to 
B = 4.414 X 10^^ G. Similar to the study of the H2 ion 5] we used a variational method with 
a simple (and unique), several parametric trial function inspired by the underlying physics of 
the problem for all range of magnetic fields studied. Our trial function can be easily analyzed 
and in contrast to other approaches our results can be easily reproduced and checked. Also 
the trial function (^ can be easily modified to explore other excited states. We showed that 
the exotic ion if|~^ exists in the states of positive z-parity Icr^, Itt^ and 16g, and does not 
exist in the states of negative 2;-parity luu, Ivr^ and ISu B = 2.35 x 10^° G — 4.414 x 10^'^ G. 
The present study complements our previous study of the ground state performed in |l| and 
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FIG. 1: Geometrical setting for the -^3^ ion placed in a magnetic field directed along the z-axis. 
The protons (marked by bullets) are situated on the z-line at distances R± from the central proton 
which is placed at the origin. 



a. 



It is evident that the state lag having no nodes is the global ground state of the system 
(pppe) (if exists) for all magnetic fields (Perron theorem). It is clear that this statement 
remains valid in general, when even the states other than studied are taken into account. 
We show for the (pppe) system in state of positive 2;-parity the electronic potential sur- 
face Et{R+, R-) always develops a minimum corresponding to the symmetric configuration 
i?+ = i?_, when the protons are situated on equal distances between each other. The do- 
main of existence of the H^'^ ion is slightly extended in comparison to the previous studies 
[l|, 1^ to be i? = 10^° — 4.414 x 10^^ G. For the case of excited states of the positive z-parity 
(l7r„, 16g) we also find a minimum in the potential surface for a similar domain of magnetic 
fields B = 2.35 x 10^° — 4.414 x 10^'^ G. A common feature for these bound states is that 
the total and binding energies grow with an increase in the magnetic field strength, while 
the internuclear equilibrium distances reduce drastically. 

For fixed magnetic field the energies of the positive z-parity states are ordered following 
the value of the magnetic quantum number m, 

E^.''' < E]7^ < e}^' , 
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as well as the equilibrium internuclear distances 

''eg ^ ''eg ^ ""eg 

This order holds also in the whole domain of magnetic fields studied (see Tables HI IIIIl and 
IIV|) . The same time for the fixed magnetic field the binding energies of the states of positive 
2;-parity are reduced slow with excitation (see Tables Hill 111 V|) . It gives us a chance to expect 
that other states of the positive 2;-parity can exist for some magnetic fields. 

For the lag state we studied the electronic potential energy surfaces for different magnetic 
fields. All those surfaces in addition to the minimum corresponding to the bound state (pppe) 
display two symmetric valleys running from the minimum to infinity corresponding to the 
"path" of the decay i^|+ -^ H:^ +P- For magnetic fields B <3x 10^^ G the /7|+ exotic ion 
is unstable towards the decay H^^ -^ H2 + p. 

Our analysis of the lowest longitudinal vibrational state and the height of the potential 
barrier (see Table ITT|l allow us to conclude that for magnetic fields B > 2.35 x 10^^ G, the 
well of the potential energy surface of the ground state lag contains at least one vibrational 
state. 

Since for magnetic fields i? > 3 x 10^'^ G the H^'^ ion is the most stable one-electron 
system, this study can be of considerable importance, in particular, in the construction of 
adequate atomic-molecular models of the neutron star atmospheres, where typical magnetic 
fields are B ^ 10^^ G or higher. A recent application of the present approach, although 
for larger magnetic fields where relativistic corrections start to be of a certain importance 



(for a discussion see e.g. W|), were used to construct a model of the atmosphere of the 



isolated neutron star 1E1207.4-5209 (see |ll|) which is based on the assumption that the 



main component of such atmosphere is nothing but the exotic molecular ion i^l"*" 12 1. 
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FIG. 2: Variational parameters of the trial function © (see text) for the state lag at the equilib- 
rium position as a function of the magnetic field strength B. Parameters ai...io are of dimension 
[a.u.]^^ and the parameters /3i,.,e are dimensionless. The parameter Ai is placed equal to 1. In the 
figure where the parameters Q8,9,io are shown the right scale corresponds to the parameter aio- 
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FIG. 3: Transition energies between the lowest three states of the H^^ ion as function of the 
magnetic field strength. 
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FIG. 4: Energy profiles along the valleys of minimal total energy for the -^3 ion as function 

of i?+,i?_, and the corresponding paths in the plane (i?-|-,i?_) for different magnetic fields: (a) 

B = 10^1 G ,{h) B = 10^2 G ,{c) B = lO^^ G , (d) B = 2.35 x lO^^ G, and (e) B = 4AU x 10^^ G. 

The position of the minimum is indicated by a bullet. The horizontal dashed line in the energy 

profile curve represents the energy of the H2 ion. 
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FIG. 4: Continuation 
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FIG. 5: Normalized electronic density distributions '^'^{x,y = 0, z)/ J '^'^(x,y,z)df and their con- 
tours for the ground state lag of the H^'^ ion along one of the valleys in a magnetic field 
B = 10000 a. u.: (a) i?-|_ = Rgq = 0.13 a. u., (b) Rj^ = 0.23 a. u., (c) i?_|- = 0.35 a. u. (near max- 
imum, i?™"^ ~ 0.36 a.u.), (d) i?-|_ = 5.0 a. u. (third proton lies outside of the figure). The positions 

of the protons are marked by bullets. Vertical axis in the figures on the left is scaled to 1 : 1000. 
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FIG. 6: Normalized electronic density distribution ^^(a:,y = 0, 2;)/J ^^(x,y, z)(ir and the corre- 
sponding contours for the first two excited states (a) l7r„ and (b) 15 g of the ion H^^ in a magnetic 
field B = 10000 a.u. The positions of the protons are marked by bullets. Vertical axis in the figures 
on the left is scaled to 1 : 1000. 
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